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Abstract 

We show that for an m-qubit quantum system, there is a ball of radius asymptotically approach- 
ing K2 _7m in Frobenius norm, centered at the identity matrix, of separable (unentangled) positive 
semidefinite matrices, for an exponent 7 = 0.5(^|-l)w .29248125 much smaller in magnitude 
than the best previously known exponent, from our earlier work, of 1/2. For normalized m-qubit 
states, we get a separable ball of radius y /3 m + 1 /(3 m + 3) x 2-( 1+ ^ m = ^3 m+1 /(3 m + 3) x 6~ m / 2 
(note that k = \/3), compared to the previous 2 x 2~ 3m / 2 . This implies that with parameters re- 
alistic for current experiments, NMR with standard pseudopure-state preparation techniques can 
access only unentangled states if 36 qubits or fewer are used (compared to 23 qubits via our earlier 
results). We also obtain an improved exponent for m-partite systems of fixed local dimension do, 
although approaching our earlier exponent as do — > oo. 

PACS numbers: 03.65.Ud,03.67.-a,03.67.Lx 
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I. INTRODUCTION AND SUMMARY OF RESULTS 



The existence of a ball of separable (that is, unentangled) multipartite quantum states 
around the normalized identity matrix, and estimates of the size of the largest such balls in 
various norms, are important for a variety of reasons. For example, lower estimates of the 
sizes of balls provide easy to compute sufficient criteria for separability of quantum states, 
as well as important tools for studying the complexity of questions about entanglement and 
multipartite quantum states. 

A series ofpapers has established the existence and provided successively better lower 
estimates S, 0, S, Q of the sizes of these balls, notably of the ball in 2-norm (Frobenius 



In this paper, we use the same general idea we used in p to obtain the best previously 
known lower estimate: the idea of considering the cone generated by tensor products of 
elements of the cone generated by a ball of separable quantum states on some multipartite 
system and elements of the cone generated by all quantum states on an additional single- 
party system. This cone will consist of separable matrices by construction; we find a lower 
bound on the radius of a ball inside it, thereby providing a lower estimate on the separable 
ball in the full system, though of smaller radius than the separable ball we started with 
on one of the subystems. By inductively or recursively combining systems in this way, we 
obtain lower estimates, dependent on the number of systems and their dimension, of the size 
of the separable ball in a multipartite quantum system. 

Here, we improve some aspects of our application of this technique, to obtain a better 
lower estimate of the size of the ball in the convex hull of the two cones (the ball-generated 
cone and the standard separable cone on different systems). When we apply the same 
inductive strategy as in p, we get a ball exponentially larger in the number of combined 
systems. For an m-partite quantum with each subsystem having dimension d , we get a ball 
of radius 



in Frobenius norm, centered at the identity matrix, of separable (unentangled) positive 
semidefinite matrices (actually we do slightly better, but with the same asymptotic expo- 
nent). For qubits (d = 2) this radius is is 

(2/3)"»/2-i > to be compared to (l/2) m ^ 1 from 
jf|. If we express it as as K2 -7 " 1 , the exponent is 7 = 0.5(|^| — 1) ~ .29248125, compared 



norm) . 



n 




(i) 
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to jj's exponent of 7 = 1/2. The non-qubit exponent is better, too, but approaches our 
earlier one as do — > 00. From this, we easily obtain a lower bound on the radius of the 
largest Frobenius-norm ball of separable normalized density matrices: for example, for m 
qubits it is (3/2) x 2~( 1+7 ) m = (3/2) x 6~ m/2 (versus our earlier 2 x 2- 3m / 2 ). A slightly 
better, but more complicated, version of our new bound lets us improve the factor 3/2 to 



3 m+1 / (3 m + 3), which rapidly approaches a/3. This gives a number of qubits below which 
NMR with standard pseudopure-state preparation techniques can access only unentangled 
states; with parameters realistic for current experiments, this is 36 qubits (compared to 23 
qubits via our earlier results). 

We also address several points not strictly necessary for obtaining these results, but which 
relate to the power and nature of our methods, and the possibilities for strengthening the 
results. Szarek [t| found the first upper bound below unity on ball size, and recently Aubrun 
and Szarek |J found an upper bound on ball size which matches (up to a logarithmic factor) 
the lower bound we obtain here for qubits, though for qudits with d > 2 there is still an 
exponential gap. One of the most natural mathematical methods for tackling this problem 
is to use a general result of F. John j^| relating the inner and outer ellipsoids of a convex 
set. We show that straightforward application of this natural method gives results weaker 
than we obtain here; weaker, in fact, than our earlier ones p. 

Our methods may appear technical; nevertheless, many of the intermediate results are 
mathematically interesting in their own right and have applications to quantum information 
problems other than the one at hand. Along the way we explain some of these, notably a 
variant proof of the result that the eigenvalues of a separable bipartite quantum state are 
majorized by those of its marginal density operators ^(|, and an example of the use of John's 
theorem to bound the radii of other inner balls of quantum information-theoretic interest, 
in this case the inner ball of the convex hull of all maximally entangled states (related to 
an application-oriented entanglement measure, the fully entangled fraction of ^j]). Many of 
our results use bounds on induced norms of various classes of maps on matrices, which we 
expect to be useful in other contexts. An appendix includes an additional bound, closely 
related to one used in the main argument, on the 2-to-oo induced norm of stochastic linear 
maps that are positive on a radius-a ball of matrices around the identity. 
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II. NOTATION AND MATHEMATICAL PRELIMINARIES 



The basic definitions and notation we use, including many elementary facts involving 

n 

cones, positive linear maps, and duality, may be found in p. Here we only review a few of 
the less-standard of these. 

We will use the term "cone" to mean a subset K of a finite-dimensional real vector space 
V closed under multiplication by positive scalars, which in addition we assume to be convex, 
pointed (it contains no nonnull subspace of V) and closed in the Euclidean metric topology. 
The dual space of a real vector space V (the space of linear functions ( "functionals" ) from V 
to R) is written V*. The dual cone to C (the set of linear functionals which are nonnegative 
on C) is C*. The adjoint of : V x -»• V 2 is <p* : V* -> V* , defined by 

(B,<P(A)) = (<P*(B),A) , (2) 

for all A 6 V\, B e V 2 *. (Here we used (B, A) to mean the value of the linear functional B 
evaluated on A.) We say a linear map <fi : V\ — > V2 is Ci-to-C 2 positive, for cones C\ C Vx, 
C2 C V2, if 4>{Ci) C C2. When C2 is a cone of positive semidefinite (PSD) Hermitian 
matrices, we will sometimes abbreviate this to "Ci-positive." 

For complex matrices M, M< denotes the transpose of the entrywise complex conjugate 
of the matrix. (The transpose itself is M f .) AoB denotes the elementwise (aka Hadamard or 
Schur) product of two matrices, defined by (AoB)ij = (A)ij(B)ij. The positive semidefinite 
(PSD) cone in the real linear space of Hermitian d x d matrices, is denoted V(d). We will 
denote by 'V" the partial order induced by this cone (X y Y iff X — Y 6 P(d))] thus 
M >l is equivalent to M G V(d). The linear space (over C) of iV x iV complex matrices is 
denoted M(N), and the linear space over the reals oi N x N complex Hermitian matrices 
is denoted 7i(N). The space of complex block matrices, K blocks by K blocks, with blocks 
in M{N), is denoted B(K,N) 

Later, we will need the following easy proposition, which follows from the fact that for 
normal (including Hermitian) matrices, A, HAH^ is the largest modulus of an eigenvalue of 
A. 

Proposition 1 Let A be Hermitian. If HAH^ < 1 then I + A y 0. 
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We use the term m-partite unnormalized density operator for a positive semidefinite 
operator 

p : H 1 ® H 2 ® H m — > H 1 ®H 2 ® ... <g> H m 
We use the term m-partite unnormalized density matrix for a matrix whose matrix elements 

p(h,i 2 , ...,i m ;ji,j2, ...,3m) 

are those of an m-partite density operator in an orthonormal basis constructed by choosing 
a fixed (ordered) orthonormal basis for each subsystem, and taking all tensor products 
® ej 2 ® • • • <g) e im of basis vectors for the subsystems. We may view this as a block 
matrix partitioned according to the value of, say, the first index; indeed, we may give it an 
m-level nested block structure (given a choice of ordering of the indices). Such a choice of 
local orthonormal bases and ordering of indices gives an isomorphism between the space of 
operators on H\ <g> ■ • • <g> H m and a space of matrices (we may occasionally implicitly identify 
these two spaces via an implicit isomorphism of this kind). 

Definition 1 Consider cones Ci C M(di), 1 < i < m. A multipartite unnormalized density 
matrix p G M(d 1 d 2 ...d m ) (corresponding to an operator 

p : H ± <g> H 2 <g> ... <g> H m — ► H X ®H 2 ® ... <g> H m ) 

is called (Ci <8> C 2 <8> ... <8> C m ) -separable if it belongs to the cone generated by the set {Ai ® 
A 2 <S> ... <S> A m : Ai G Cj, 1 < i < m}. We call this the separable cone, S(Ci, C 2 , C rn ). 

This is trivially equivalent to the recursive definition: 5(Ci,C2, ...,C m ) is the cone gen- 
erated by the pairs A <g> B with A G S(C 2 , C m _i), B G C m , and S(Ci) := C 1 . 

When Ci (for 1 < % < m) are the PSD cones V(di), (Ci <8> C 2 <8> ... <8> C m )-separability is 
the standard notion of separability of multiparty unnormalized density matrices. 

We will use various norms on spaces of matrices or operators, including the Frobenius or 
2-norm \ \A\\ 2 := V tr A^ A, the 1-norm ||A||i := tr vALA, and the operator norm ||^4||oo '■ = 
maxilla,! | =1 In the definition of the operator norm, we used vector norms (written as 

1 1 ■ 1 1) on the input and output spaces, which we will take to be the Euclidean norms induced 
by our chosen inner products on these spaces. In general for linear operators : V — > W 
and norms || • ]| T and || • ]| w on V, W respectively, we will write 

11011™ := max \\(f>(x)\\ u ; (3) 

xeV:||x|| T =l 
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this is the operator norm induced by the choices r, u for norms on V, W. Also, when 
(f) : M(K) — > M(N) is Hermitian preserving, we write 4> H for 0's restriction to Hermitian 
matrices (i.e. to have domain H(K) and range Ti(N)). These details are motivated by the 
fact that key technical results of our paper involve the relationship between norms (induced 
by various choices of matrix norms on the input and output matrix spaces) of linear maps 
: M(K) — > M(N), and similar norms of 4> H . 

Finally, a note on our usual choices for naming dimensions, which should help make things 
clearer below. When considering a multipartite Hilbert space Hi ® H 2 <8> ■ • • <E> H m , we use 
d\, d 2 , d m for the dimensions of Hi, H m , and d for the overall dimension n™-^. When 
we consider combining a ball cone and a PSD cone (as described in the introduction and in 
more detail below), we let the ball cone be in a space of d 2 x d 2 Hermitian matrices, and the 
PSD cone in a space of d\ x d\ Hermitian matrices. When we consider linear maps between 
matrix spaces, we usually use the somewhat unnatural choice that M(d 2 ) (or H(d 2 )) is the 
input space, and M(di) (H(di)) the output space. When we consider an m-partite system 
where all the subsystems have the same dimension, we use do for the dimension of a local 
system and d for the total dimension d™. 



ILL MAIN RESULTS 



We begin with some key definitions; then we give an outline of the proof of our main 
results, followed by the detailed proof. 

Definition 2 If X is a bipartite density matrix viewed as an element of B(di,d 2 ), so that 
its blocks X l) i are in M{d 2 ), and if : M(d 2 ) — ► M(di) is a linear operator then we define 

1 ^(X 1 ' 1 ) 0(X X ' 2 ) ... ^(X 1 ^) N 
^(X 2 ' 1 ) 0(X 2 - 2 ) ... 0(X 2A ) 

^(X^' 1 ) 0(X dl < 2 ) ... 0(X dl ' dl ) y 
A simple result characterizing separability, but one fundamental to our argument, is: 

Lemma 1 Suppose that the cone C(d 2 ) C H(d 2 ) C M(d 2 ) . Then X is V(d 1 ) ® C(d 2 )- 
separable iff <f>(X) >z (i.e. is positive semidefinite) for all stochastic C{d 2 ) -positive linear 
operators : M(d 2 ) — > M(d 1 ). 
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For the proof, see (f| . 

With these, we can sketch the proof of our main result, which applies to a tensor product 
of systems of dimensions d\, d 2 , d n . It is a recursion relation for a radius a n such that 
all matrices within (or at) Frobenius norm distance a n of the identity are separable (i.e. 
P{d\) ® P{d 2 ) ■ ■ • P(c? n )-separable): 



a n < a 



n-l< 



d„ 



(5) 



l ^2(i-aL 1 /(nr = -M))K-i) + i ' 

Proof- outline: 

1.) Begin by letting d 2 in Lemma ^ be the total dimension U^~^di for our set of systems 
and C(d 2 ) be the separable (i.e. P{d\) ® ■ ■ ■ <E> P(<i n _i)-separable) cone for these systems, 
and d\ of the lemma correspond to d n for our n systems, so the lemma says X is separable 
if and only if: 



1 (piX 1 * 1 ) <P(X l > 2 ) ... <P(X 1 ^) ^ 
<p(X 2 > 1 ) 0(X 2 ' 2 ) ... 0(X 2 ^) 



y o 



(6) 



v 0(X dl!l ) 0(X dl - 2 ) ... 0(X dl < dl ) J 
when (f)(S(d 2 , d n )) C P{d\) and 0(7) = I. 

2.) Since the ball Ball(a n _i) of radius a n _i around the identity is separable by hypothesis, 
the set of stochastic operators that are positive on that ball is no smaller than those positive 
on the separable matrices, so X is separable if © holds for all such 0. Let X = I + Y , Y 
Hermitian and traceless; by Proposition ^ X is separable if 



IIWIIoo<l. 



(7) 



3.) For stochastic <p with 0(Ball(a n _i)) >z 0, we easily show ^(M)^ < (l/a n _i)||M|| 2 
when M is Hermitian, while for M traceless but not necessarily Hermitian we obtain 



IMAOHoo < A = (l/a n _i) y/2(l - a 2 jd x d 2 ■ ■ ■ d n _!) I \M\ \ 2 . 



(8) 
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4.) We bound the LHS of (J7J) with elementary norm inequalities (for typographic clarity, 
inside the norm delimiters, we omit the curved braces that otherwise delimit block matrices): 



0(Y 2 ^ 



0(Y 2 > 2 ) 



< 



<p(Y dul ) (j){Y d ^ 2 ) 

Mr^Woo H0(^ 2 )Hoo . 
^(y^lU ||0(r 2 ' 2 )||oo • 

H^Y^IU ||0(F dl ' 2 )iioc • 

dll^lla A||y^|| 2 ... 
Alalia a-UWY^W, ... 



X\\Y d ^\ 



\\\Y d ^\ 



</>(Y 



H0(^ Ml )llc 
||0(F 2 ^)||c 



M Y d ^)\\ 

A||F 1A || 2 
M\Y 2 ^\\, 



a n-l\\ I 



< 



(9) 



where we used the bounds from step 3.), along with the fact that Y's offdiagonal blocks may 
be made traceless by local transformations without affecting its separability or entanglement, 
in the last inequality. 

5.) We prove an upper bound on ||0_b||2^oo for maps 4>b '■ Z \— > B o Z, and evaluate it 
in the case that I?'s matrix elements are equal to a constant on the diagonal, and another 
constant off the diagonal. Calling this upper bound /ig, we have ||</>b(^)||oo < A*b||-^||2) we 
apply it to the last expression in step 4.) to get: 



°n-lll^ /rl,1 ||2 A||y 



1,2 I 



'I4i 



2.1 1 



-1 \W2,2\ 



Ally 1 - -M2 

A||y 2 ' di || 2 



Ally 



A||y di ' 2 || 2 ... a-^Hy^'^i 



< 



yx, 



i.i i 



yx, 



1.2 I 



ly^ib i|y 2 ' 2 ii 2 



|y Mi ll 2 
|y 2 ' di ii 2 



||ydi,i|| 2 ||F^ 2 || 2 ... ||Y*.*|| 2 

= A*Bi|y||2 



(10) 



s 



By step 2, then, /xell^lh < 1 5 i- e - ll^lb < implies separability of X = I + Y. We 
have 



B 



a n-l ^ • • • ^ 

A a n _! ... A 



y X A ... a n _ x y 
and we will show that ps works out to be 



/ A 2 K-l) + a 2 
V d n 

Using the expression we will derive for A gives that 

ll*1|2< 



d r , 



A 2 (d n -1) + 1 



d r , 



(11) 



(12) 



(13) 



\| 2(1-^/^-^X^-1) + ! 

guarantees separability of / + Y, establishing (jSJ). 

We will apply our results also to balls of separable normalized states, using the following 
result taken over from J(| (where it is Proposition 7). This proposition is derived using 
"scaling," i.e., considering all ways of writing a matrix p as a positive scalar times the sum 
of the identity and a Hermitian perturbation, and minimizing the 2-norm of the perturbation. 



Proposition 2 Define p{p) as the maximum of ||A||2 over all A such that there exists an 
a > for which p = a(I + A). Let p be a normalized (tx p — 1) density matrix. Then the 
following three statements are equivalent: 

1. p(p) < a. 

2. trp 2 < l/(d-a 2 ). 

3. \\p-I/d\\ 2 < a/^d{d-a?). 

Corollary 1 Let a be a lower bound on the size of the m-partite separable ball around the 
identity matrix, d be the dimension of the m-partite Hilbert space. If an m-partite normalized 
(i.e. unit trace) density matrix p : Hi®- ■ -®H m — > ■ -®H m satisifes \\p—I/d\\2 < ^, 

where d = dim[H\ ® ■ ■ ■ ® H m ), then it is separable. 
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(The proposition actually gives the (negligibly) tighter statement with yj d(d — a 2 ) in the 
denominator.) 

We now embark on a more detailed presentation and proof of our results, beginning with 
some definitions. 

Definition 3 Let G(N, a) C H(N) C M(N) be the cone generated by hermitian N x N 
matrices of the form {I + A : 1 1 A| I2 := (tr(AA^)^ < a}. 

Let : M(d 2 ) — > M(d\) be stochastic. Consider the maximum "contraction or dilation 
ratio" of on Hermitian operators, 

max MA)\U\\A\\ 2 . (14) 

Hermitian A 

Note that this is equal to max[| J 4.[[ 2 =i,A Hermitian I \4>{A) | |oo, and therefore equal to 

110*112-00. (15) 

Definition 4 Define j(di, d 2 , a) as the maximum, over stochastic maps : M(d 2 ) — > M(di) 
that are positive on G(d 2 ,a), of 1 1 0^ 1 1 2—^00 ■ 

Note that we used here, not itself. 

Proposition 3 Let Hi,H 2 have dimensions di,d 2 . If an unnormalized density matrix p : 
Hi ® H 2 — > Hi (g) H 2 satisfies the inequality \\p — I\\ 2 < l/ r y(di,d 2 ,a) then it is V(di) ® 
G(d 2 , a)-separable. 

Proof: Let p = I + A, A Hermitian; by Lemma 1, we are looking for a bound on ||A|| 2 
that ensures, for any stochastic G(d 2 , a)-positive linear operator (i.e. 4>{X) >z for all 
X e G(d 2 ,a)), that 4>(I + A) h 0. 0(7) = /, so 0( I + A) = I + 0(A); H^A)^ < 1 
will ensure 0(7 + A) >z (cf. Proposition [TJ). Since ||0(A)| l^/l |A|| 2 < -y(di,d 2 ,a) from the 
definition of 7(^1, d 2 , a), ||A|| 2 < l/ r y(di,d 2 ,a) ensures this. ■ 
In order to make good use of this proposition, we need a bound on the value of 7(^1, d 2 ,a). 
Proposition 0] below, together with Proposition EJs bound on the parameter A that appears 
in Proposition HJ provides it. Obtaining this bound on _7 is the technical heart of our results, 
and the improvement in this bound over that found in p is the source of the better exponent 
in the lower bound on the size of the separable ball we obtain in the present paper. We 
begin with a definition and an easy lemma. 
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Definition 5 Define X(di,d 2 ,a) as the maximum, over all stochastic maps <p '■ M(d 2 ) — > 
M(d\), positive on G(d 2 ,a), and over all traceless X e M(g? 2 ), of H^pf)!^/] |X|| 2 . 

Lemma 2 If <p : M(d 2 ) — > M(d\) is a stochastic G(d2, a) -positive linear map with < 
a < I, and 0(7) = I e M{d\), then 

< « _1 H^ll2 for allX e H{d 2 ). 

Proof: G(d 2 , a)-positivity of a stochastic means | |0(A)| ^ < 1 for all Hermitian | |A| | with 
||A|| 2 < a; since H^A)!^ is homogeneous in ||A|| 2 , it will achieve its maximum on such A 
where ||A|| 2 = a, implying | |$(A)| U/l | A| | 2 < 1/a. ■ 
We now proceed to our key bound, on 7(6^, d 2 , a). 

Proposition 4 Suppose a > l/d 2 . Then 



a 2 X 2 (d u d 2 ,a)(d 1 - 1) + 1 
7(rfi,rf 2 ,a)<a W . (16) 

Proof: Let A E B(di,d 2 ) be a Hermitian d\ x d\ matrix of d 2 x d 2 blocks Call the 

2- norms of the blocks afj, and the operator norms of the blocks afj, and define A 2 and 
A°° as the matrices with these elements. Similarly, call the matrices whose elements are 
||0(A( JJ '))||{ 2)OO }, ${2,00}^ (y^ e p rom i se no t to square any matrices named A or so this 

notation is unambiguous.) 

Note that ||A 2 || 2 = ||A|| 2 . Also, note that 

MA)\\ao< Halloo, (17) 

by an elementary norm inequality (the operator norm of a block matrix is bounded above by 
the operator norm of the matrix whose elements are the operator norms of the blocks of the 
original matrix). $oo is a matrix with nonnegative entries. Its diagonal entries ^(^'^Hoo 
are bounded above by a _1 ||A( v )|| 2 by Proposition^! which applies since the diagonal blocks 
of A are Hermitian. The offdiagonal blocks are not in general Hermitian, but they may 
be made traceless via a unitary "local transformation" (acting only on the index specifying 
which block) which has no effect on the matrix's separability or entanglement. This is 
because one of its (unnormalized) "reduced density matrices," is the matrix of traces of its 
blocks, and the reduced matrix may be diagonalized by a local transformation. 
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So for the offdiagonal entries i ^ j we have = | | ^ < A(g?i, tZ a , a) | | b 
by the definition of A. In other words, using < for the ordering in which A < B means B — A 
is (entrywise) nonnegative, we have 

$°° < a~ x L o A (2) , (18) 

where L is the matrix with l's on the diagonal and a\(d,2, di,a) =: rj in all offdiagonal places. 
Therefore (since the operator norm is monotonic in the ordering <), the maximal contrac- 
tion/ dilation ratio on Hermitian matrices, i.e. the 2-to-oo induced norm (on Hermitian 
matrices), of the completely positive map A a : X 1— > a~ 1 L o X taking M(K) — > M(K) is an 
upper bound on H^Ua-nxi- The induced norm of A a is a _1 ||A||2_ >0 o, where A : X \— > L o X; 
we evaluate it via the following Lemma. 

Lemma 3 Let <pB be the linear map from TC(n) to Ti{n) defined by <pB '■ X 1— > B o X, for 
some Hermitian B. Then 

||0s||2^oo= max y l Cy , (19) 



where C is the n x n matrix with elements Cij = |Sy 
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The lemma states that the 2-norm-to-oo-norm induced norm of the positive map defined 
by the Schur (elementwise) product with B for some fixed Hermitian B, is just the maximum 
value of a quadratic form over a simplex, the matrix of the quadratic form being the one 
whose elements are the absolute squares of B's. This lemma has independent interest; we 
defer its proof and a discussion of other applications to Section IIVI 

Recall the abbreviation aA(d 2 , d±, a) =: 77, and note that the premise of the Proposition 
we are proving implies 77 > 1. We have 

l|A||Loo= max (l-tftEtf+^O) 2 

= (i-v 2 )T,y" + v\ (20) 

i 

where we used J2i Hi = !• Since rj > 1, this is maximized where J2i vf is minimized, i.e. with 
each yi = l/d\. The maximal value is (1 — r] 2 )/di + rj 2 , and thus 
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Since (as argued before Lemma EJ) |A 2 | | 2 ^oo is an upper bound on | \<p H | |2-kx>, this 
gives the desired result. ■ 
Remark: The ease with which we were able to use Lemma |3] in the above proof was due 
to the simple form of the matrix L which took the role of B. The problem of maximizing 
a general quadratic form with nonnegative matrix, over the simplex, is NP-hard as one can 
reduce Max-Clique to it (this is apparently well-known, cf. [3] or jl^]). 

To make further use of this in evaluating 7(^1, o^, a )i we need an estimate for A(g?2, d%, a). 
The following proposition provides one. 

Proposition 5 



A(g? 2 , di, a) < -y 2(1 - ^) . (22) 
ay «2 



This plays the role that Proposition 6 did in [6J, but while that proposition did not as- 
sume (J) stochastic, and established that for all <fi whose 2-to-oo-induced norm on Hermitian 
operators is at most 1, the induced norm on all operators is at most a/2, the present propo- 
sition adds the assumption of stochasticity, and computes the maximum induced norm for 
the class of stochastic G(d2, a)-positive maps acting on traceless matrices, rather than all 
matrices. In fact, using Proposition 6 of p for the bound on A and the rest of the argument 
as in the present paper, we could have obtained the same exponent in our bound on ball 
size as a function of number of systems m. 

Proof: We need good bounds on the 2-to-oo induced norms of G(d2, a)-positive maps 
: M(d 2 ) — > M{d\j. Since it will turn out that these do not depend on d\, we will use 
d in place of g?2 throughout the discussion. We consider normalized matrices in G(d,a), 
which are expressible as p = I/d + A for some traceless Hermitian perturbation A, and 
recall from Proposition |2] that these are precisely those normalized p for which ||A||2 = 



\\p— I/d\\2 < a/yd(d — a 2 ). G(d, a)-positivity is equivalent to positivity on these normalized 
matrices (since they generate the cone G(d, a) by positive scalar multiplication). The latter 
is equivalent to the condition 



(j)(I/d + A) y whenever 1 1 A| \ 2 < a/yjd(d - a 2 ) . (23) 
Using Proposition 1, for stochastic 4> this is equivalent to 



lAHoo < l/d whenever ||A|| 2 < a/^d(d-a 2 ) . (24) 
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For Hermitian traceless A, | |0(A) | \oo/\ | A| I2 is homogeneous of degree zero in A, and therefore 



||0(A)|U/||A|| 2 < (l/d)/(a/y/d(d-a?)) = 



(l/a)^l-a 2 /d . (25) 



To extend this to arbitrary, not necessarily Hermitian, traceless matrices B write B in terms 
of traceless Hermitian and traceless antiHermitian parts as B = X + iY . Then 

Halloo < I loo +11*001100 



< a-Vl-a 2 M||X|| 2 +||F|| 2 ) 



< aT 1 y/l - a 2 /dV2b , (26) 
where the second inequality is ()25|) and the last is elementary Euclidean geometry. Thus 



II^IU/H^lb < aT l p{\ - o?/d) . (27) 

■ 

Incorporating the upper bound of Proposition explicitly into Proposition 0] gives 
Proposition 6 



7(d x , d 2 , a) := max max [^(A)!!^/! |A|| 2 

<p Hermitian A 



a 



_! 2(l-a 2 /rf 2 )(d 1 -l) + l 



Using this bound in Proposition El gives: 



(28) 



Proposition 7 Let Hi,H 2 have dimensions di,d 2 . If an unnormalized density matrix p : 
H\ ® H 2 — > Hi ® H 2 satisfies the inequality 



ll "- f|l '^ a V 2(l-»'/W-D + l <29) 
then it is V(di) ® G(d 2 , a) -separable. 

We may apply this proposition inductively or recursively, in various ways, to obtain 
bounds on multipartite separability. In the following the induction proceeds as in by 
tensoring one additional PSD cone P(d n ) with a cone generated by a ball of separable states 
in P(di) <S> ■ ■ ■ <E> P(d n _i), of radius a n _i, obtained in the previous inductive step. The 
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induction begins with the base case of a bipartite separable ball of radius one (in 2-norm) 
around the identity (from 5]). From Proposition [7| we have the recursion relation: 



a n < a n-lf 



(In 



(30) 



l \\2(l-al_ 1 /(TVt- 1 %))(d n -l) + l • 

This allows for easy numerical calculation of a n . When we have a total of m systems each 
of dimension do, we have: 



0-n < Ovi-1, 



dn 



\2(l-a 2 n _ 1 /d n - 1 )(d -l) + l 



For qubits, this is 



0"n < 0>n-l\ 



(31) 



(32) 



3 - a2_i/2"_ 1 " 

Using the weaker bound X(di,d 2 ,a) < a~ l \/2 from [6( gives a weaker but easily solved 
recursion relation: 

Proposition 8 Let Hi,H 2 have dimensions di,d 2 . If an unnormalized density matrix p : 
H\ (gi B.2 — > Hi <S> H% satisfies the inequality 



\\p - J|| 2 < a x 



di 



2(di - 1) + 1 



(33) 



then it is V(d\) <g> G(d2, a) -separable. 

This gives a worse bound, but asymptotically the same exponent for the number of systems: 



Corollary 2 If an m-partite unnormalized density matrix p:Hi®---® H m — > Hi ® • • • ® 

H m satisfies 

Wr-'W^y^) (34) 

then it is separable. 

While for large do, Corollary |21 is asymptotically the same as the bound (1/2)-/- from Q, 
for qubits it gives the notably better 

Corollary 3 If an m-qubit unnormalized density matrix p : Hi® - ■ -®H m — > Hi® - ■ -®H m 

satisifes 

\\p-I\W < (2/3) m/2 - 1 (35) 

then it is separable. 
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In fact, we may explicitly solve the recursion (|31|) exactly, obtaining: 
Theorem 1 If an m-qudit unnormalized density matrix satisfies the inequality 



1 lp ~ J| ' 2 ^ r " := V (2d-l)"-£p-l) + l • (36) 



then it is separable. 
For qubits, we have: 



in— 1 



Qn+1 

_ 2 - 7 m ( 37 ) 



3 n + 3 



with 7 = 0.5(g| - 1) « .29248125, compared to 7 = 1/2 from [6j. 

In an earlier version of the present paper we obtained the same exponent, but a slightly 
worse overall expression, because we did not exploit local transformations to render the 
offdiagonal blocks of A in the proof of Proposition 0] traceless, and so had to use a slightly 
worse contraction bound X'(di,d,2,a) that applies to all matrices, not just traceless ones. 
Since this bound A' may prove useful in other situations, we include it and its proof in an 
appendix. Subsequently, Roland Hildebrand 14] obtained the same asymptotic exponent 
but a slightly larger ball for the m-qubit case, via an argument exploiting the fact, special 
to the case of m qubits, that the local cones are already ball-generated (aka Lorentz) cones. 
In the proof above, we exploited the ability to render the offdiagonal blocks of A traceless 
by local transformations, improving the bound to agree with Hildebrand's in the qubit case, 
but also improving it for the case of m (io-dimensional systems (and indeed, in general). 

Although Corollary El gives a ball with k = 3/2, as mentioned above, we see that the 
present paper's improved bound on X(d\, d%, a), as embodied in (|3*T|) gives a larger ball, with 
a prefactor «;(m) asymptotically approaching a/3. For tripartite separability of unnormalized 
states, Proposition [7| gives a ball of radius \/4/5 around the identity (a result also noted by 



Hildebrand), larger than our previous result of y 1/2. 

Using Corollary El and Proposition |21 for m qubits, we obtain a lower bound on the radius 
of the largest normalized separable ball of 2" m (2/3) m / 2 - 1 , i.e. (3/2) x Q- m ' 2 = (3/2) x T~^ m 



with rj = (l/2)(ln6/ln2) ~ 1.292481. Using the stronger recursion we get J^rWZ ( 1+ T m ) = 



3"+3 

course of investigating the volume of the separable states relative to 
all normalized states, Szarek (|7] Appendix H) obtained a lower bound of Q~ m / 2 on the 
radius of a related, but larger "symmetrized" set S, the convex hull of S U —S). In general 
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case such symmetrization can substantially increase the inner radius. Indeed, in the case of 
the (i-dimensional simplex, the inner ball has radius of order 1/d compared to 1/ \fd for its 
symmetrization (which is the unit sphere in Zi-norm). Szarek also obtained the first upper 
bound below o(2~ m ) on the radius of balls inside the normalized separable m-qubit states: 
it is o{2~ rrm ) with the exponent rj = 1 + (1/8) log 2 27/16 w 1.094361). Recently Aubrun and 



Szarek 



improved this, obtaining an upper bound for for the symmetrized set of separable 



normalized states of qubits (which contains the separable states) of: 



C 0V /mlogm6- m/2 . (38) 

r-i 

The constant Co is equal to VSC\, where C\ (which appears in a crucial lemma of |8j) can 
be chosen to be 1.67263, and can probably be chosen smaller. The asymptotic exponent 
for this expression matches that in our lower bound for the case of qubits, though with the 
logarithmic pref actor. On the other hand for d > 3 the inner radii in the unsymmetrized 
and symmetrized cases are of different order. Indeed, it is easy to prove that that the 
unnormalized separable radii r(d\, ...,dk) > r(Di, Dk) if di < D^l < i < k. In 8j, 
Aubrun and Szarek also state an upper bound of (d (d + l))~ m / 2 (up to a similar prefactor) 
for the normalized symmetrized qudit case, corresponding to order ((d+ l)/d)~' m l 2 for the 
unnormalized ball around /. This should be compared to our results for the unnormalized 
ball which are of order ((2d - l)/d )~ m/2 - While both of these give (3/2)" m / 2 in the case of 
qubits, the Aubrun-Szarek exponent (with a constant base such as 2) approaches zero as do 
grows, while ours does not (approaching, instead, —1/2). Thus in the case of do > 3 there is 
still an gap between our result and their upper bound, and it is an interesting open problem 
to close this gap. Notice that it had been proved in p] that the radius of the maximum ball 
inside the normalized real-separable m-qubit states is 0(2~ m ) = 0(l/d) (indeed, it is exactly 



1/ Jd(d — 1) = 0(1/ d) for general real-separable multipartite states). We also showed in 



that the bipartite _b.e states have in-radius 1/^-1) (resolving a question raised, 
for example, in [15|, where the d = 4 case was proved). The 0(1/ d) results correspond 
to a ball of radius order unity of unnormalized real-separable or bipartite separable states, 
compared with one that (from Szarek's upper bound) must shrink as an inverse of a power of 
dimension in the general unnormalized multipartite case. This provides another example of 
a dramatic difference in the behavior of entanglement in the bipartite versus the multipartite 
situation. 



17 



IV. A 2-TO-oo INDUCED NORM BOUND FOR SCHUR PRODUCT MAPS 



In this section, we prove Lemma|2l which was used in proving Proposition0]in Section lTTTl 
It states that the maximum 2-norm-to-oo-norm contraction/dilation ratio for the positive 
map defined by the Schur (elementwise) product with B for some fixed Hermitian B, is just 
the maximum value of a quadratic form over a simplex, the matrix of the quadratic form 
being the one whose elements are the absolute squares of S's. 

Lemma 3 Let <f>B be the linear map from TC(n) to Tt(n) defined by 4>b '■ X i— > B o X , for 
some Hermitian B. Then 

||0s||2^oo= max y*Cy , (39) 



where C is the n x n matrix with elements Cij = \B{j 



|2 



Proof: To show this, we use the basic fact (see e.g. |16|) that for vector spaces (finite- 
dimensional, for simplicity) V, W equipped with norms 1 1 • | |v, 1 1 • | \w, and using the notation 
II " \ \v*, || ■ \ \w* f° r the norms dual to || • ||y, || • \\w, for a fixed linear map T : V 

muL\\T(X)\\ w /\\X\\ v = max ||r*(y)|| v V||F||^ . (40) 

Using the facts that 0b is its own dual {(f) b = <Pb)^ the 2-norm is its own dual norm, 
and the operator norm is dual to the 1-norm, we obtain that max||x|| a =i ||0-BpOI|oo = 
max||y|| 1=1 ||0 B (y)|| 2 . We proceed to evaluate the latter. 

Since ||«/>bC^)||2 is increasing in ||^||i, the maximization can be extended to the convex 
set {Y : \\Y\\i < 1}, and since ||0b(F)||2 is convex the maximum will occur at an extremal 
point of that set. The extremal points of the ball ofnxn Hermitian matrices with 1-norm 
at most 1 are the rank-one projectors (pure states) X whose matrix elements are XiX*, for 
some normalized (J^iX? = 1) vector x. For such X, 

(<j> B (X))i3 = B » X i X *j ■ ( 41 ) 



Hence 



II0bWII2 = EI 5 -I 2 I^I 4 + E^AI^I 2 I^I 2 

= l^| 2 N 4 + E l%l 2 N 2 N 2 • (42) 
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Defining y as the vector in R+ with yi := \xi\ 2 and the matrix C by CV,- = l-By] 2 , this 
expression is just y t Cy, and we are to maximize it over y e R+ such that Y^iVi — 1; 
establishing the lemma. ■ 
Digression: Completely positive maps of the form considered in Lemma E3 are useful in a 
variety of contexts in quantum information theory. Simplest, perhaps, is their appearance 
in the most general representation of "partial decoherence" processes in some basis. The 
relevant mathematical fact here is that the set of completely positive maps T such that there 
exists an orthonormal basis for which T(ejeJ) = e^ef, or equivalently all states diagonal 
in that basis are fixed points of the map, is precisely the set of maps X i— > B o X with B 
Hermitian and having ones on the diagonal. These maps are doubly stochastic, implying 
that the output density matrix is "more disordered" than the input density matrix, meaning 
its eigenvalues are ^ by those of the ioput oW y ^ (0, Theorem 7.1). 

Another application is an alternative proof of a fact due to Nielsen and Kempe [lo| . 
that the vector of decreasingly ordered eigenvalues of a separable bipartite mixed state is 
majorized by that of either of its marginals (reduced states): "separable states are more 
disordered globally than locally." The proof uses the well-known fact, useful in a variety 
of contexts both within and outside of quantum information, that the (necessarily PSD) 
matrices AA^ and A^A have the same eigenvalues. Equivalently, a quantum state (even an 
unnormalized one) 

R = Y, «V f (43) 

i 

has the same eigenvalues as the Gram matrix of the (not necessarily normalized!) vectors 
v l (the matrix whose ij element is the inner product (v l ,vi) = v l 'tP), as one sees by letting 
A in the above fact be the matrix whose i, k element is the k-th coordinate of v 1 in some 
orthonormal basis. A separable state R (even unnormalized) has a representation of the 
form (}4*3|) with v 1 = x % ® y l , where we may take \\y l \\ = 1 without loss of generality. Its 
eigenvalues are therefore those of the Gram matrix G with elements 

G tj = ^VW . (44) 

The marginal state on the first factor is J2iX l x l \ whose eigenvalues are those of H whose 
elements are 

Hij = x i] x j . (45) 
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But 

G = B o H , (46) 

where B is the Hermitian PSD matrix, with ones on the diagonal, whose elements are 
Bij = y^yi . Therefore (by the abovementioned fact that the eigenvalues of the output of 
a doubly stochastic map applied to a Hermitian operator are majorized by those of the 
Hermitian input), G's eigenvalues are majorized by if's, proving the statement. ■ 

V. COMPARISON WITH AN APPROACH VIA JOHN'S THEOREM 

A celebrated result of Fritz John |9j is a natural tool for approaching this problem, so we 
verify here that our methods provide stronger results than one can get by straightforward 
application of John's theorem. John's theorem gives a shrinking factor such that, when the 
smallest ellipsoid covering a convex set is shrunk by that factor, it fits inside the set. This 
is interesting in itself; and if we know the ellipsoid, then we can obtain (from its shortest 
axis) a ball that fits inside the set as well. 

A. The inner and outer ellipsoids, the coefficient of symmetry, and John's theorem 

Let S be a closed compact convex set (of nonzero measure, i.e. generating the vector 
space V) in a real vector space V of dimension D. Let E ont be the least- volume ellipsoid 
containing S. Let Sintered be S translated so that the center of its E out is at the origin. 
Define the "coefficient of symmetry" of S as the largest "shrinking factor" < a < 1 such 
that for every x in Sintered, — otx is also in Sintered- John's result states that if we shrink 
the least- volume covering ellipsoid, -Eout, by multiplying it by a factor yoc/D the resulting 
shrunken ellipsoid is contained in S. Note that when a set S is symmetric under the action 
of a compact group G, so are E in (S) and E out (S). 

B. Application of John's theorem to the set of normalized separable states 

Every ellipsoid in the normalized quantum states is a set of the form: Lq := {p : Q(p — 
I/d, p—I/d) < 1}, for a quadratic form Q that is strictly positive semi definite on the positive 
semidefinite matrice, and for every such form Q, Lq is an ellipsoid. 
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Proposition 9 Let E = {p : Q m i n (p — I /d, p — I /d) < 1} be the minimum-volume ellipsoid 
covering the m-partite separable normalized density matrices. Let d = n^di, where d{ are 
the local dimensions. Then 

{P : Q mi „(p - IMP- I/d) < = ^§= C S . (47) 

Proof: We first calculate the coefficient of symmetry, by noting that the quantification over 
x in the definition of a can be restricted to x extreme in Sintered '■= S — I/d, i.e. shifted 
versions 7r — I/d of pure separable states tt. Let x be an arbitrary extremal state; we find 
the largest a such that —ax G Centered- That is, we seek the largest a such that 

— a(ii — I/d) G S — I/d, i.e. 

(1 + a)I/d-an G S . (48) 

The LHS of (|4*H|) (which we'll call A a ) has unit trace for all a, and is PSD (certainly a 
necessary condition for its separability) as long as a < l/{d — 1). With this value of a, it 
becomes: 

A-^V-). (49) 

We now show that A is separable, so the coefficient of symmetry is l/(d — 1). Since it is 
separable, it is equal to x^x^ ® x^x^ <S> • ■ ■ <S> x^x^ 71 ^ for some normalized vectors 
x ( m ) g ff m p or e ach p g {i ; ; m | i e t {^r }ie{i,...,d p } be a complete orthonormal basis with 
first member x&\ Then, since i a ... t m x h x il ® ' ' " ® ^iT 4™ = becomes 

(ii, 12, ...,i m )^(l,l, ...,1) 

This expresses A as a convex combination of separable pure states, demonstrating A's sep- 
arability. 

Since a = l/(d — 1) and D = d 2 (the dimension of the real linear space TC(d) of d x d 



Hermitian matrices) we have Ja/D = (l/d)Jl/(d — 1). John's theorem then gives (|47|). 



Remark: The smallest ball B covering S is centered at I/d and has radius y(d— l)/d. 
This follows from the easy fact that the pure separable states (indeed all pure states) lie on 
the boundary of this ball, which by unitary invariance therefore contains all the normalized 
states, including the separable ones. If this ball were E then (|4*7|) would give us a ball of ra- 
dius 0(d~ 3 / 2 ) inside the Hermitian matrices. When the system consists of m do- dimensional 
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systems, this is 0(d,Q^ 3 ^ 2 ^ m ). For qubits, this would have the same exponent as the results in 
p, though it would still be less good for general d Q (where |(| gives 2~( 1//2 ) m <io m ). The results 
we obtain elsewhere in this paper always have a better exponent, though it converges to the 
exponent of our earlier result as dimension grows. E is in fact not a ball (we thank Stanislaw 
Szarek for pointing this out to us). Still, the above result establishes that straightforward 
application of John's theorem does not give us better results than or the techniques we 
use in the other sections of the present paper. The largest ball we can straightforwardly 
get via John's theorem is the largest ball in the shrunken minimum-volume ellipsoid, whose 



radius is l/(dy/d — 1) times the length of the least principal axis of the covering ellipsoid E. 



This must be no larger than d 3 / 2 = (l/(dy/d — l))J(d — l)/d, for if the least principal axis 



of E were larger than the radius ^J(d — l)/d of the smallest covering ball B then E could 
not be minimum-volume. 

We note that a natural approach to obtaining E itself is to use some of the more ele- 
mentary aspects of the methods exposed in noting that E out =: E (and E- m ) must be 
invariant under the action of conjugation by local unitaries JJ\ ® Ui ® ■ • • ® U m , E must be 
a ball when restricted to each irrep of this action; finding the radii of each of these balls 
determines E. 

We also note that in the bipartite case, our maximum ball Bailed) of the radius = 



l/yd(d — 1) in the Frobenius norm is, in fact, also the maximum- volume ellipsoid inscribed 
in S. Indeed, it had been proved in |5| that this ball Bailed) belongs to the convex compact 
set of normalized separable bipartite states; on the other hand it is easy to show that Bailed) 
is the maximum-volume ellipsoid inside the (larger) convex compact set of all normalized 
bipartite states. 

Remark. Group symmetry can easily be used to compute the coefficient of symmetry for 
other convex hulls of orbits of interest in quantum information theory (and thus when E out 
can be computed, one gets lower estimates of the inner ball's radius via John's theorem). 
For example: 

Proposition 10 Let F denote the convex hull of all normalized "maximally entangled 
states" of a bipartite system with local dimensions n (overall dimension d = n 2 ), i.e. the 
convex hull of the orbit of the state tc := \l/\Ed G B{n, n), where 

* = (1/Vn)$>i® ei , (51) 
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under the action of U(n) x U(n) on B(n, n) as conjugation by local unitaries: (u, v) G 
U{n) x U{n) acts as: X \-> (U(u) <g> V(v ))X{U{uft <8> V{vY), U,V being standard n x n 
matrix representations. The coefficient of symmetry of F is l/(d—l). 

Proof: As before, by symmetry it suffices to find, for a single extremal state tt 6 F (for 
which we choose 7r as defined in the Proposition), the largest a such that (|38|l holds, with 
F substituted for S. Exactly as before, we get a < l/(d— 1) = l/(n 2 — 1) necessary for 
positivity. We want to show that when a = l/{n 2 — f), the state 

(f + a)I/d -aix = (l/(n 2 - 1))(J - n) =: R (52) 

is not only positive but in F. To this end we use the Choi/Jamiolkowski isomorphism, 
and view the matrix R as associated with a map T. I is the Choi matrix of the map 
Z : X i— * (tr X)I, i.e. the projector onto the one- dimensional subspace of matrices spanned 
by the identity, while 7r is the Choi matrix of 1/n times the identity map id : X i— > X. 
Therefore R is the Choi matrix of 

' -(Z--id) (53) 



n 2 — 1 v n 

For every finite dimension n, there exists at least one orthogonal basis U^i G 0, ...,n 2 — f 
for M(n) with Uq '■— I, and all Ui unitary. (For example consider the basis {P k S l : k,l £ 
{0, ...,n — 1}}, with P the diagonal matrix whose j-th diagonal element is a/ -1 for some 
primitive n-th root of unity u, and S is the matrix with elements Sy = %+i) mo d «,j mod n\ 
the gene ral question of which such bases exist is considered in ^|.) It is easily verified (cf. 
e.g. [lflj ) that for any such basis the map Z may be written 

n 2 -l 

Z : X h+ (1/n) J] • (54) 

i=0 

Therefore, with the notation T4 for the map X 1— > AXA\ R is the Choi matrix of 

j n 2 -l 

7 2 u E T ^ • (55) 
Since T4 has Choi matrix n(I ® A)7r(J ® A'), (|55|) implies 

1 n 2 -l 

i? = - T -— ^ (/ ® ^)tt(J ® , (56) 

which expresses it as a convex combination of local unitary transforms of tt, as desired. ■ 
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VI. APPLICATION TO THERMAL NMR STATES AND PSEUDOPURE STATES 



In many interesting experimental or theoretical situations, the system is in a "pseudopure 
state": a mixture of the uniform density matrix with some pure state ir: 

p £i7r := err + (1 - e)I/d , (57) 

where d = di,..d m is the total dimension of the system. For example, consider nuclear 
magnetic resonance (NMR) quantum information-processing (QIP), where d = 2 (the Hilbert 
space of a nuclear spin), and m is the number of spins addressed in the molecule being used. 
As discussed in more detail below, the initialization procedures standard in most NMRQIP 
implementations prepare pseudopure states. 

Using Corollary Q with b a lower bound on the unnormalized 2-norm ball around I, p e ^ 
is separable if 



d-1 



e<(b/d)^j—^<b/^d(d-l), (58) 
For m c?o-dimensional systems (so d = do" 1 ), this implies the (negligibly loosened) bound 

e < b/d m . (59) 



e < — ath , (60) 



Since we have established in this paper a bound of b = (d Q /(2d - l)) m/2 -\ we obtain 

1 

d™ /z+L (2d~- l)™/ 2 " 1 

This is an exponential improvement over the result in j^J (the qubit case is in P]) of e < 
1/(1 + d 2 ," 1-1 ), and indeed over our results in p, although as d — > oo the improvement in 
the exponent of m over that in [6] goes to zero. 

In liquid-state NMR at high temperature T, the sample is placed in a high DC magnetic 
field of strength B. Each spin is in a thermal mixed state, with probabilities for its two 
states (aligned (f) or anti-aligned (j) with the field) proportional to e ±/3/iB , where (3 = 1/kT 
with k Boltzmann's constant, /i the magnetic moment of the nuclear spin. For realistic 
high-T liquid NMR values of T = 300 Kelvin, 5 = 11 Tesla, rj := PfiB w 3.746 x 10~ 5 < I. 
Since e ±v m 1 ± T], the probabilities are m (1 — rj)/2, p[ ~ (1 + f])/2. Thus the thermal 
density matrix is approximately 

1+2 Q 
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(with each qubit expressed in the |t)?|l) basis). The highest-probability pure state 
of independent distinguishable nuclear spins, has all m spins up and probability about 
(1 + r]) m /2 m ps (1 + mT])/2 m . Standard pseudopure-state preparation creates a mixture 

(l-e)//2- + e|T---T)(T---T|, (62) 

where 

e = V m/2 m . (63) 

of this most probable pure state and the maximally mixed state, by applying a randomly 
chosen unitary from the group of unitaries fixing the all-spins-aligned state. With rj ps 
3.746 x 1CT 5 , this implies that below 36 qubits, NMR pseudopure states are all separable, 



compared to the ~ 23 qubits obtained in 



and the ~ 13 qubits one gets from the bound 
in [3j]. Since we have not shown that the bounds herein are tight, with our assumed rj 
even at 36 qubits there is no guarantee one can prepare an entangled pseudopure state by 
randomization. We remind the reader, also, that if such a state existed, there would still 
be no wayof partitioning the qubits so that the state exhibited bipartite entanglement; as 



noted in |6], the results of [5] imply that for the parameters used above, one needs m = l/i] 
qubits (about 26, 700 for our rj) before the pseudopure state obtained from the thermal state 
by the randomization procedure described above fails to satisfy pf's sufficient criteria for 
bipartite separability with respect to any partition of the qubi ts into two sets. 

Schulman and Vazirani's algorithmic cooling protocol [20( shows that it is, in theory, 
possible to prepare any entangled state of n qubits from polynomially many (in n) thermal 
NMR qubits, although the overhead is discouraging. The question of just how many qubits 
are required by means possibly simpler than algorithmic cooling is also of interest. One can 
gain some information about this using our results, by applying Corollary ^ to the initial 
thermal density matrix of an NMR system. For the initial thermal density matrix p of 
flBTjl .we have: 

I \P ~ I/d\\l = { -—^ « W/2 m • (64) 

This should be compared to the separability condition obtained by using the relation (|3T)|) 
for a n , and Corolary^ Numerical comparison shows that 17 qubits are required before this 
bound is exceeded (rather than the 36 required for the pseudopure state prepared from this 

25 



thermal state). (Our earlier bound allowed only the weaker statement that for fewer than 
14 qubits, no entanglement exists in the thermal state p.) 
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APPENDIX A: CONTRACTION BOUND FOR STOCHASTIC BALL-POSITIVE 
MAPS ON ALL MATRICES 

In this section, we state and prove a contraction bound from 2-norm to oo-norm (i.e. a 
bound on the induced operator norm) for stochastic, ball-positive maps on all matrices. It 
is slightly more involved to prove than the one for maps on traceless matrices used in the 
body of the paper, but although we ultimately did not need it for the present paper, we 
present it here in the hope that it may find uses elsewhere in quantum information theory 
or mathematics. 

Definition 6 Define A'(g?i, d,2, a) as the maximum, over all stochastic maps <fi : M(d,2) — > 
M(di), positive on G(d,2,a), arid over all X e M(d 2 ) , of \\(p(X)\\ 00 /\\X\\2. 

Proposition 11 



X'(d2, di,a) = 




Proof: Recall from (J25|) that for Hermitian traceless 

\\<j>(A)\U\\A\\ 2 < (l/a)y/l-a?/d . (A2) 
To extend this to arbitrary, not necessarily Hermitian traceless, matrices consider: 

M = c(I/Vd) + B , (A3) 
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with B traceless but not necessarily Hermitian. To bound ||0(M)| loo/I |M|| 2 it suffices by 
homogeneity to bound it for ||M||2 = 1, i.e. defining \ \B\\i =: b, for c 2 + b 2 = 1. Writing B 
in terms of Hermitian and antiHermitian parts as B = X + iY , we have: 

IWAOiu < c/Vd+||0(x)|| oo + ||$(y)|| oo 



< c/^ + a-Vl-a'Mll^lh + \\Y\\*) 



< c/\fd + ar l Jl-a?/d^/2b 



(A4) 



where the second inequality is by (jA2|) and the last is elementary Euclidean geometry. 
Defining 



7 := a~V2(l - a 2 /d) 



(A5) 



we maximize the RHS of (|A4|) over c, 6 such that c 2 + fe 2 = 1 (i.e. | \M\ [2— 1). We obtain 



l + 7 2 (i' 



7 2 d 



l + 7 2 d 



2^ ' 



and hence a maximal value for the RHS of 



Substituting our definition for 7 gives 



||0(^OI|oo < \/2/a 2 — l/d . 



(A6) 



(A7) 



(A8) 



Thus an upper bound on ||0(M)| loo/I |M|| 2 for arbitrary M and G(d,2, a)-positive stochastic 



<p (which is to say on \(d 2 , di, a)) is y2/a 2 — l/d 2 

For the lower bound portion of the proposition , we exhibit a G(d, a)-positive stochastic 
map map r for which | \t(X) \ \oo/\\X\ [2 = \J^z - |- We begin by defining a family of stochas- 
tic maps parametrized by \x > 0, acting on Hermitian matrices. For N > 4, d\ > 2 we define 
r by specifying r(I) = I, and: 

\ 



( 1/V2 
-1/V2 











V 



/ 





-1 



fJ>cr z 
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/o 

-1/V2 

l/y/2 

\- ■ ■ ■ 



( 



1 ... 



1 ... 



J 



Dots indicate the matrices are to be filled out with zeros. I and the two input matrices 
given above are mutually orthogonal in trace inner product; on the orthocomplement of 
their span, r is taken to map everything to zero. Call the input matrices above Z and X (so 
that Tfj,(Z) = a^,r M (X) = a x ). r M extends to antiHermitian matrices homogeneously, due 
to its Hermiticity preserving property, so that r^{iZ) = ifia z , r^(iX) = ifia x . (The names 
a x , a z are chosen for the output matrices because the usual Pauli matrices that go by these 
names appear in the upper left-hand 2x2 blocks of our a x ,a z , and are padded out with 
zeros.) 

For Hermitian traceless B, the maximal value of 1 1^(5)1 |oo/||-S||2 will occur where B = 
cZ + bX. Then t^(B) = fi(ca z + ba x ) = c 2 + b 2 )a a , where a a is some matrix which 
has a 2 x 2 Hermitian upper left diagonal block with eigenvalues ±1, and is zero elsewhere. 
■■ fiVc 2 + b 2 , and since ||S|| 2 = v 'c 2 + b 2 , | |t(.B) | |oo/| | 2 



Hence \\r(B)\\ a 
for to be G(d, a)-positive requires that 



My)|| M /||y|| 2 < a -yi-a 2 /d 



hold for all traceless Hermitian Y (cf. ()25|) ). so we must have 



ji. However, 



(A9) 



V» ^ a 1 y 1 — a 2 /d . 



(A10) 



We choose \i equal to the RHS here; then the inequality (jA9|) holds for all Hermitian Y, as 
required for G(d, a)-positivity. 

Now, we consider the not-necessarily-traceless matrix Y = -^=1 + (P/\/2) ( X + iZ) . Then 
\\Y\\2 = V a ' 2 + P' 2 , which we set equal to one WLOG. Now, 



T 



(Y)\\oc = a/Vd+ (/3/V2)||0(X) +i</>(Z)\\ O0 
= a/Vd + (3fi\\a x + ^0-21100 



a/Vd + (3a~ l \ll 



l /dV2. 



(All) 
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The last equality uses just the definition of // and the result 1 1 cr^. + zcr^loo = 2. The latter is 
easily obtained by noting that 

/ 



i 1 
1 -i 




(A12) 



1/V2 

This vector [i\/2, a/2]* has Euclidean norm 2, which is therefore a lower bound on the 
operator norm of a x + ia y ; since the Frobenius norm upper-bounds the operator norm, and 
is equal to 2 in this case, the operator norm is 2. Define 



i := a-y2(l - a 2 jd) = ^2/a 2 - 2/d . (A13) 

Then we have 

\\r(Y)\\ 00 = a/Vd + (3i (A14) 



and the same argument used to obtain (|A7|) as the maximum of (|27)j) yields J 1/d + r )" 1 as 



the maximum here. Substituting the definition of 7' gives a maximum of y2/a 2 — 1/d for 
A, which matches the previously obtained upper bound. ■ 
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